Finite size effects on pion spectra in relativistic heavy-ion collisions 
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We compute the pion inclusive transverse momentum distribution assuming thermal equilibrium 
together with transverse flow and accounting for finite size effects and energy loss at the time of 
decoupling. We compare to data on mid-rapidity pions produced in central collisions in RHIC at 
yijvjv ~ 200 GeV. We find that a finite size for the system of emitting particles results in a power- 
like fall-off of the spectra that follows the data up to larger pt values, as compared to a simple 
thermal model. 



PACS numbers: 25.75.-q 



I. INTRODUCTION 



In recent years, strong evidence has been found in favor 
of the production of matter composed of the fundamen- 
tal QCD degrees of freedom in central Au-|-Au collisions 
at the relativistic heavy- ion collider (RHIC) Some 
of the quantitative understanding about the signals sup- 
porting these findings are associated to intermediate or 
high Pt phenomena in particle spectra. Among these, the 
depletion in the production of large pt charged hadrons 
and neutral pions, as compared to extrapolations of data 
on p+p collisions to RHIC energies Q , suggests that the 
fragmenting partons that give rise to final hadrons, loose 
energy on their way through the dense partonic matter 
believed to be formed during the early stages of the reac- 
tion. This phenomenon has been named jet quenching Q 
and has received stronger experimental support from the 
results of the d-f-Au mode at RHIC 5]. Moreover, the be- 
havior of the proton to pion ratio for intermediate pt 2 
GeV 6| suggests that an important mechanism of hadron 
production at RHIC is the thermal recombination of free 
quarks JJ during the evolution of the collision. 

In spite of the success of this conceptual framework to 
interpret RHIC data, there is an ingredient commonly 
overlooked: the fact that all these phenomena take place 
during small time scales and consequently within small 
volumes. 

Although not commonly considered, it is certainly true 
that small size effects are important in the description 
of statistical systems. For example, finite size effects 
have influence on the late-stage growth of nucleated bub- 
bles during a first order phase transition Also, in 
the statistical hadronization model, the states used for 
the description are in general different to the observed 
asymptotic ones but such difference is not an issue un- 
less the volume of the system is small, typically of order 
0(10 im)^ 9]. Finite size effects are also known to influ- 
ence the interpretation of the correlation lengths in Han- 
bury Brown- Twiss anal ysis in the context of relativistic 
heavy- ion collisions uM ■ 

Finite size effects on particle spectra have also been re- 



cently investigated 0, 0, 0| , where these were referred 
to as boundary effects. The scenario considered in the last 
two of these studies resorts to a mean field treatment of 
the interactions affecting the most abundantly produced 
particles in the reactions, namely pions. It was argued 
that prior to kinetic freeze-out, the pion mean free path 
is smaller than the pion interaction range and therefore 
that during this stage, pions are still strongly interacting, 
mainly attractively, among themselves. The system re- 
sembles more a liquid than a gas, with a surface tension 
acting as a reflecting boundary. When assuming ther- 
malization at freeze-out, the pion spectrum can be com- 
puted from the sum of the square of the wave function 
in momentum space weighed by the thermal occupation 
number for each state. 

Qualitatively, the behavior of particle spectra includ- 
ing this finite size effects deviates from a simple exponen- 
tial fall-off at high momentum since from the Heisenberg 
uncertainty principle, the more localized the states are in 
coordinate space, the wider their spread will be in mo- 
mentum space. In terms of the discrete set of energy 
states describing the pion system, this behavior can be 
understood as arising from a higher density of states at 
large energy as compared to a calculation without bound- 
ary. 

In this work we compute the transverse momentum 
distribution for pions assuming thermal equilibrium to- 
gether with transverse flow and accounting for finite size 
effects at decoupling. By comparing to data on pion spec- 
tra on Au-I- Au coUisions at ^/snn = 200 GeV (Qt il5|j , we 
show that for temperatures and collective transverse flow 
within the commonly accepted values, the transverse mo- 
mentum distributions can be described up to larger val- 
ues of Pt as compared to a simple exponential fall-off. 

The paper is organized as follows: In Sec.^ we present 
the basics of the model to compute the pion spectra in- 
cluding finite size effects. In Sec. IIIII we compute the 
thermal pion spectrum including the effects of transverse 
flow together with the description of finite size effects. 
We carry out a systematic study of the behavior of the 
spectra when varying the parameters involved. In Sec. II VI 
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FIG. 1; Invariant pion distribution as a function of pt for 
different combinations of the parameters T and (3 for a fixed 
value R = 8 fm. The normalization is arbitrary. 



FIG. 2: Pion rapidity distribution for T = 120 MeV, (3 = 0.6 
and R = S fm. The normalization is arbitrary. 



we compare the calculation to data on pion transverse 
distributions from Au+Au collisions at y'SNN = 200 
GeV and show how this model pushes the thermal com- 
ponent of the spectrum to larger pt values as compared to 
a simple exponential description. Finally we summarize 
and conclude in Sec. |3 



II. THE MODEL 

When the system of particles can be considered as con- 
fined and their wave functions as satisfying a given con- 
dition at the boundary just before freeze out, the en- 
ergy states form a discrete set. The shape of the volume 
within the confining boundary is certainly an important 
issue in case the model aims to describe transverse mo- 
mentum particle distribution at all rapidities. However, 
in this work we will restrict ourselves to describe spec- 
tra at central rapidities. We thus consider a scenario 
in which the system of particles of a given species is in 
thermal equilibrium and is confined within a sphere of 
radius R (fireball) as viewed from the center of mass of 
the colliding nuclei at the time of decoupling. This time 
needs not be the same over the entire reaction volume. 
Nevertheless, in the spirit of the fireball model we con- 
sider that decoupling takes place over a constant time 
surface in space-time. This assumption should be essen- 
tially correct if the freeze-out interval is short compared 
to the system's life time. 

The solutions that incorporate the effects of a finite 
size system have been found in Ref. 13]. They are given 



as the stationary solutions to the equation 

subject to the condition 

ib{\r\=R,t)=0, (2) 

and also finite at the origin. The normalized stationary 
states are 

^ , ,(r t) ^ 1 / 1 y/^ 

R Jl+3/2{knlR) \rEnl J 

i1r„'(r)J,+i/2(fcnir)e-'^"'* (3) 

where Ji, is a Bessel function of the first kind and Y/m' is 
a spherical harmonic. The parameters kni are related to 
the energy eigenvalues Eni by 

Eli = kli+m' (4) 

and are given as the solutions to 

Jl+l/2{knlR) = . (5) 

Thus, the eigenfunctions for the bound system of parti- 
cles are given explicitly by Eq. ©. These are analytical, 
albeit given in terms of Bessel functions. 

III. THERMAL SPECTRUM 

The contribution to the thermal particle invariant dis- 
tribution from a state with quantum numbers {n, I, m'} 
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FIG. 3: Invariant pion distribution as a function of pt for two 
values of the parameter R and fixed values T = 120 MeV and 
P — 0.6. The normalization is arbitrary. 
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FIG. 4: Nuclear suppression factor Raa for neutral pious from 
Ref. parametrized in the interval from 1 to 10 GeV. The 
parametrization is chosen as to become a constant for large 
Pt values. 



is given by 



E- 



da 
(2^ 



\k,ii ■u)f{kni ■ v)Wnim'ip,r), (6) 



where Wnim' (p, r) is the Wigner transform and /(/c„i • v) 
the thermal occupation factor of the state tpnim' , respec- 
tively. The four- vectors v'^ and represent the collective 
flow four-velocity and a four-vector of magnitude one, 
normal to the freeze-out hypersurface a, respectively. 

For the system of bound bosons, the Wigner transform 
of a given state is defined as 



r' d r' 

In order to consider a situation where freeze-out happens 
at a fixed time and within a spherical volume of radius 
i?, the unit four-vector can be chosen as u*^ — (1, 0). 

To keep matters simple, we consider a thermal occu- 
pation factor of the Maxwell-Boltzmann kind, namely 



flkrd • u) = e 



(8) 



where T is the system's temperature. The four-vector 
is parametrized as 

v^=-^{l,^r), (9) 
and we choose a radial profile for the vector v such as 



v = /3 



(10) 



where the parameter [3 represents the surface expansion 
velocity. Notice that, the radial profile for the expan- 
sion velocity considered in Eq. H10|l can be taken as the 



transverse component of the fiow only for the description 
of central rapidity data. Correspondingly, the gamma 
factor is given by 
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(11) 



Nonetheless, in order to continue to keep matters as 
simple as possible, and be able to analytically perform 
the integrations in Eq. ((BJ, we will instead consider that 
the gamma factor is a constant evaluated at the average 
transverse expansion velocity, namely. 
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VI - (3/3/4)2 
We take the four-vector fc^^ given by 

(£'„;, k„i). 



(12) 



(13) 



and choose k„; || p. This choice is motivated from the 
continuum, boundless limit, where the relativistically in- 
variant exponent in the thermal occupation factor be- 
comes ^{E — p • v). 

Gathering the elements described in Eqs. to H13|l . 
it is straightforward to evaluate the integral in Eq. 10. 
Summing over all the states, the invariant thermal dis- 
tribution is given by 
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FIG. 5: Invariant pion distribution as a function of pt for 
R — 8 fm, T — 120 MeV and (3 = 0.6 compared to data on mid 
rapidity positive pions from central collisions at ^sjvjv = 200 
GeV measured by PHENIX. The normalization of the model 
theoretical curve has been fixed from a fit to the data. 
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FIG. 6: Invariant pion distribution as a function of pt for 
T = 120 MeV, /3 = 0.6 for values of J? = 8 fm (upper curve) 
to i? = 14 fm (lower curve) , compared to data on mid rapidity 
neutral pions from central collisions at -yijvjv = 200 GeV 
measured by PHENIX. The normalization of the curves has 
been taken from the fit to the tt^ spectrum. 



where the factor {21 + 1) comes from the degeneracy of 
a state with a given angular momentum eigenvalue I and 
A/" is a normalization constant. 

Figure ^ shows the invariant pion distribution as a 
function of pt for different combinations of the param- 
eters T and (3 and a fixed value i? = 8 fm. Figure QJi 
shows the distribution for two different temperatures 
T = 100, 120 MeV and (3 = 0.5. As expected, the effect 
of increasing the temperature is to increase the inverse 
slope of the distribution. Figure^ shows the distribution 
for two different values of the surface expansion velocity, 
/3 = 0.5, 0.8 and T = 100 MeV. The effect of increas- 
ing the expansion velocity is also to increase the inverse 
slope of the distribution. Figure |2 shows an example of 
the kind of rapidity distributions obtained in the model. 
Notice that the distribution is not as broad as data seem 
to indicate (see for example Ref. ^^■) This is a limi- 
tation of the spherical symmetry assumed in the model 
emphasizing the need of restricting its applicability to 
the description of central rapidity data. 

Figure 13 shows the distribution for two different values 
of the confining radius R = 8, 10 fm for fixed values of 
T = 120 MeV and (3 = 0.6. The effect of increasing the 
confining radius is to increase the value of pt for which 
the distribution has a point of inflexion towards a concave 
shape. This is to be expected since in the limit where the 
confining volume goes to infinity, the value of pt for this 
point of inflexion should also go to infinity as the system 
becomes boundless. 



IV. PION SPECTRA 

In order to test the model, we compare the theoretical 
distribution to data on mid-rapidity pions produced in 
RHIC central collisions at ^/^^ = 200 GeV [1 E|. 
Rather than performing an overall fit to find the optimum 
parameters, we fix them to values within the commonly 
accepted ones to describe freeze-out conditions. 

To include the effects of parton energy loss, we resort 
to parametrize the data on Raa for neutral pions as mea- 
sured by PHENIX We use the the expression 



and find from a fit to the data on Raa the values 
A = 0.18, B = 3.17 and C = 0.20. This parametriza- 
tion is taken in such a way that for large pt the nuclear 
suppression factor remains a constant. This is shown in 
Fig.H 

First we describe data for the low pt part of the spectra. 
This is done in Fig. |31 where we show the invariant pion 
distribution as a function of pt for R ~ 8 fm, T = 120 
MeV and /3 = 0.6, including the energy loss effect, com- 
pared to data on 7r+ from PHENIX ^J. The normaliza- 
tion of the model theoretical curve has been calculated 
from a fit to the data. We notice from Fig. O that the 
curve does a very good job describing the data for all 
values of Pt. In contrast, a calculation where no effects 
of a finite size are included, and thus the wave function 
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of a given state is simply a plane wave, does not describe 
the data over the considered range when we use the same 
values for T and (3 as for the case of the calculation with 
a confining volume. This is also shown in Fig. |S1 where wc 
notice the more steep fall-off of the curve without finite 
size effects. 

Next we describe data in the range 1 GeV < pt < 10 
GeV. This is done in Fig. where we show the invari- 
ant pion distribution as a function of pt for T = 120 
MeV, f3 = 0.6 and i? = 8 — 14 fm, including the en- 
ergy loss effect, compared to data on neutral pions from 
PHENIX . The normalization of the curves has been 
taken from the fit to the tt"*" spectrum. The model does 
a good description of data for values of the parameter 
R between 8 and 10 fm up to ~ 3 GeV from where 
the curves start deviating from the data that shows a 
steeper fall-off. This failure is however expected since 
for large pt the leading particle production mechanism is 
the fragmentation of fast moving partons, some of which 
fragment outside the fireball region and thus are not in- 
fiuenced by the confining boundary that the rest of the 
particles experience within the fireball. 

V. SUMMARY AND CONCLUSIONS 

In this work we have computed the pion transverse mo- 
mentum distribution in a thermal model including trans- 
verse expansion and accounting for the fact that the sys- 
tem of particles is produced within a small volume up to 
freeze-out, in the context of relativistic heavy-ion colli- 
sions. The finite size has been introduced by requiring 
that the system's wave functions vanish at the bound- 
ary but remain finite at the origin. The net effect is a 
broadening of the particle transverse momentum distri- 
bution with respect to a simple thermal model without 
a boundary. The physical origin of this broadening is 
the Heisenberg uncertainty principle since as the states 
are more localized in coordinate space, their spread is 
wider in momentum space. Energy loss has also been 
considered by a parametrization of the experimentally 
measured nuclear suppression factor Raa where we re- 
quire that this becomes constant for large pt values. 



The model does a very good job describing data for low 
to intermediate pt in the range < pt ^ 3 GeV for com- 
monly accepted values of the temperature T = 120 MeV 
and the surface transverse expansion velocity P — 0.6 
(corresponding to an average transverse expansion veloc- 
ity (/3) ~ 0.5) when the system's radius R at freeze-out 
is in the range 8 fm < i? < 10 fm. This is best seen when 
comparing the model to data on charged pions produced 
in central Au -I- Au collisions. The description of data for 
Pt in the range 3 < < 10 is not as good as can be seen 
from a comparison of the model to the spectrum of neu- 
tral pions produced in central Au -t- Au collisions. This 
behavior is not surprising though, since for large pt it is 
well known that the leading particle production mecha- 
nism is the fragmentation of fast moving partons, some 
of which fragment outside the fireball region and thus 
are not influenced by the confining boundary, or mean 
field, that the rest of the particles experience within the 
fireball. 

We thus conclude that the thermal features of the spec- 
trum, which are obtained from the low pt region, are best 
highlighted when accounting for the finite size of the sys- 
tem at decoupling. It is of course interesting to check 
the consistency of this treatment when extended to the 
case of fermions, and in particular protons. Finally, it is 
worth to note that finite size effects should also play a 
role during the partonic phase of the collision jlTI] and in 
particular that they can also influence the recombination 
scenario. This kind of analyses is for the future. 
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